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Abstract: The asymptotic iteration method (AIM) is an iterative technique used to find exact 
and approximate solutions to second-order linear differential equations. In this work, we employed 
AIM to solve systems of two first-order linear differential equations. The termination criteria of 
AIM will be re-examined and the whole theory is re-worked in order to fit this new application. As 
a result of our investigation, an interesting connection between the solution of linear systems and 
the solution of Riccati equations is established. Further, new classes of exactly solvable systems of 
linear differential equations with variable coefficients are obtained. The method discussed allow to 
0^ ' construct many solvable classes through a simple procedure. 
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I. INTRODUCTION 

There are a limited number of exactly solvable systems of differential equations that are known in the literature [H-[l]. 
' For example, the system of two constant-coefficient differential equations: 

0'2 =C(/)i +d(j)2 (1) 

has a general solution depends on the roots of the characteristic equation 

\^ -{a + d)\ + {ad-hc) = Q. (2) 

■ The purpose of the present work is two-fold: (I) to introduce a new technique based on the asymptotic iteration 
' method 3-0 to analyze the exact solutions of linear systems 

0| (/)2 =wo(a;)0i -I- /3o(a;)02, (3) 

ON . 

, where Ao(a;), sq{x), ujq{x) and po{x) are continuous and differentiable functions, and (2) to add new classes of exactly 
solvable systems to the known ones [Jj-f^. 

k>l \ The asymptotic iteration method (AIM) is an iterative technique introduced originally Q to find the exact and 
j_j ■ approximate solutions of second-order linear homogeneous differential equations Q. Since then, the method has 
d ' enjoyed a number of interesting applications, particularly, in relativistic and non-relativistic quantum mechanics 

[T]| . In their study of Dirac equation, Hakan et al Q modified AIM to study Dirac equation by analyzing a system 
of linear differential equations, however, their work was limited to this particular application. In the present work, we 
extend AIM to deal with different systems of first-order linear differential equations in general, and to give a detail 
structure useful in generating exactly solvable systems. 

The paper is organized as follows: in next section, we re-examine AIM to deal with a system of linear differential 
equations. In section III, we used AIM to provide a simple recipe to find the exact solution of a system of two 
constant-coefficient differential equations ([T]). In section IV, we give a necessary condition for polynomial solutions 
of linear systems of differential equations ([3]). In section V, we establish an interesting connection between the exact 
solution of linear system ([3]) and the exact solution of a Riccati equation; namely, for every exactly solvable Riccati 
equation, there exist at least one system of linear differential equations that is exactly solvable. In section VI, some 
applications and tables of results are reported. The conclusion is given in section VII. 
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II. ITERATIVE SOLUTION OF LINEAR SYSTEMS 



We consider a system given by ([3]), direct differentiation with respect to x yields 



62 =U}l4>l+ Pl4>2 



where 



Ai = Ao + Ao + so^o 
si = So + AqSo + so/Oo 
cji = + iOoXa + paUjQ 
Pi = Po + ^oso + Po 

In general, the n"* and {n + 1)*'* derivative of (jlj are: 

^(ri+1) ^ ^^^^ _|_ ^^^^ 

and 

01 = A„+i0i + S„+i(/)2 

02 = + Pn+102 

respectively, where A„_|.i, s„+i, d^n+i and are computed recursively by: 

An+l = '^n + -^riAo + SnCjQ, 
Sn+1 = + A„So + S„Po, 
W„+l = Cj'j + CJnAo + Pnl^Oi 
Pn+1 = P'i + t^nSo + PnPO- 

The ratio (/)i"+^V0i"^^^ then reads 

^, ^ A„+i(0i + |^02) 

rfa; ^ A„(0i + f^02) ■ 

If, for sufficiently large n > 0, we have 

then ([9]) reads 

and consequently, we have 

= Ci exp ^^^^i) = CiXn exp (awo + Ao) dt 
where we used ^ and (fTU]) . Further, using ([S]), we have 

01 + q;02 = Ci exp {auo + Ao) 
which can be solve for 0i to yields, using ([3]), the first-order differential equation of 02 

02 + (ujoa - po) 02 = ujqCi exp I / {uJoa + Ao) I . 



The exact solution of this first-order differential equation is easily obtained; namely 

(^2 = exp (po - woa) dt^ Ci+Cx J (^cjq exp (^J {Xq- po + 2woa) dr^ ^ dt 
where Ci and C2 are the constant of integration. Consequently, we have 

^1 = Ci exp (awo + Aq) dt^ - a(f>2- 

Equivalently, if we consider instead the ratio (^^^^"^^ l'^^^^^^ ^'^'^ for sufficiently large n > 0, we have 

Wn+l _ _ ^ 



Pn+1 Pn 



that 



and 



51 =exp 



j (Ao - So/3) dx^ C2 + C[ J (^so exp (^J {po-Xo+ 2so/3) dx^ ^ dx 



<p2 = -/3<j)i + C[ exp {/3sq + po) dx 

where again C[ and C2 are constant. The above discussion establish the following theorem: 
Theorem 1 The general solution of the system: 

(f)[ =Xo{x)(t)l + So{x)(t)2 



is given by: 

4>2{x) =exp (^J {po{t) - aujoitj) dt 
<pi{x) =Ci exp {au)o{t) + Xo{t)) dt ] - a4>2{x 

if for sufficiently large n > 



i)'2 =UJo{x)<j)i + Po{x)(j)2 

C2 + C1 j LOoit) exp 



(Ao(r) - po(t) + 2a;o(T)a) dr 



a 



Sn+l 



Xn Xn 



+ 1 



Equivalently, if for some n > 0, 



then the general solution is given by: 



1^ _ '^n _ 

Pn Pn+1 



Mx)=e^p(^J {Xo{t)-pso{t))dt^ C2 + C[J (^so(i)exp(^y" (po(r) - Ao(r) + 2so(t)/3) dr 

(j)2{x)=-(3Mx)+C[exp (^J iPsoit) + po{t)) 



where 



)dt 



An+i = A^ + A„Ao + s„a;o, 

Sn+l = s'n + A„So + S„Po, 



: + UJnXa + PniOo, 



Pn+1 = Pn+ ^nSO + PnPO- 
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TABLE I: Exactly solvable systems iji^ = Ao(a::)<j!>i + so{x)(l>2, and <j!>2 = a;o(a;)<^i + po{x)(p2, for Xo{x),so{x),u![}{x) and po{x) 
arbitrary functions. Here, a, b, c, and d are constants and n — 1,2, ... . 



Xo{x) 



so{x) 



LUo{x) 



Solution 



bc-(n-l)'' + (n-l)d 
(d~n + l)x 



h — C'l 



ad-{a + d)(n-l) + (n-iy^ 



e(l + d-n) 1 d4 



+ C2S" 



/>! =Cl 
/>2 = Cl 



1 + Q — 



2+a+d-2ri -1 



2+a+d-2n J 



Cs. 



(d-n + l)(a-n + l) 



Cl 
Cl 



„ , , . „ \X ■■ + C2[ ^]x"~ 

(l + a-n)(l + tj-n) ]_l + a + d-n , i^^^ri-l 
b(2 + a + d-2n) ^^-f- 



(l+g-n)^ 



bc+(n — l)a — (n — 1)~ 
(a — n + l)^; 



/>i =Ci 
A2 = Cl 



a^ + !)c-2(n-l)a + {n-l)^ 
c(l + a — n) 



+ C2 



a^ + i)c-2(n-l)a + {n-l)^ 



+ C2S; 



n — a — 1 -I 
n-1 



Note that ([22|) and (|23p can be written in more convenient way as 



and 



T — T— On+l — A„4-iS„ — An-Sn+l — U, 



— A„+i — tJ„+ip„ — UJnPn+i — U. 



Pn+1 Pn 

Example 1: Consider the following linear system of differential equations 



(26) 



(27) 



a 




b 




X 




a; 


h 


c 








X 


h H 


X 


h 



(28) 



where a, 6, c, and d are arbitrary constants with ad~ cb ^ 0. It is straightforward, using (j26p . to show that 

= '[[{■m'^ -mia + d) + ad-bc), n=l,2,... (29) 

m=0 



X 



.)2 -{n-l){a + d)+ad-bc = Ooi a - ^^-("-1) +("-1)'^ 
_ ^ ("-!) g^]-^(j consequently the system, for n = 1, 2, . . . 



Clearly, (5„ = 0, if {n — ly — {n — l){a + d) + ad — be = or a — d-n+i ■ Furthermore, we have, using ((22|) . 



6c- (n - 1)2 + (n - l)d 



{d — n + l)x 

d . 



x 



X x 



(30) 



has the general solution 



--Cl 



be 



bc+il + d- n)2 

c{l + d- n) 
6c+ (1 + rf - n)2 



n ~ d — 1 



x'^+T+d^ + C2X 



ri-1 



(31) 



for n — 1,2, ... . Similar cases are reported in Table 1. Further examples of exactly solvable systems are reported in 
Table II. 
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TABLE II: Exactly solvable systems ([T| for different Xo{x),so{x),ljo{x) and po{x) using Theorem 1 where a,h,c,d, and ai,i = 
1, . . . , 4 are arbitrary real numbers. Here n — 1,2, . . . and Ea[z] is the exponential integral function Ea[z] = t''"'e~^^dt. 



\o{x) 



so{x) 



uo[x) 



po{x) 



Solution 



a + 



^1 cd+ax ^""^ J {cd + axY^ 



'^dx"-'Jjlj:^dx]-C2Cx" 



{cd + ax)x"-\CidJ T^^dx + C2) 



a + 



b-n + l 

d X 



-^V 14. b — n + ax ' V I J J + b — n + ax)^ ^ ' ^ d 

h = -x"-\l + b-n + ax){CidJ ^f^l^-^dx + C2). 



a + 



(b-n + l) 



Ci [ nl. "I > + (1 + fe - n)x"-i f ,"7 " w rfa:) + CC2X'' 

-^L(l + £i — n + ax) ^ -' J (1 + 6 — n + ax}^ 

""-1(1 + & - n + aa;)[^iii±^ 



1 + 

+ i> — n + ax) 



^dx + C2] 



+ ■ 



02 



as + 



14 + - 



4i = Cia: "4 ^ 



[-(^ + ai)x\ + a4a:£; ^4+°.^ [-(^ + ai)x]) + C2a;"-\ 



III. SOLUTION OF LINEAR CONSTANT-COEFFICIENT SYSTEMS 



Theorem 2 The constant- coefficient first-order linear system 

4>'i =\o4>l + Sq4>2 

where Ao,so,t^o aJ^c? po are constants, has the general solution 



4>i ^Ci 



ujQa 



UJQ 



^(Ao+woaja; _|_ ^^g(po-woa)a; 



where a is given by: 



So 



2a;n 



or a 



(po - Ao) - V(Ao - PoY + 4woSo 



(32) 



(33) 



(34) 



Proof: From the AIM sequence ([25]), we have for the constant-coefficient linear system that 

An+l = A„Ao + SnLUo, 
Sn+1 = A„So + SnPO, 
= i^n^O + PntOo, 
Pn+1 = ^^nSO + PnPO- 



So + apo 



Consequently, the ratio ([22)1 now reads 



which yield the quadratic equation 



Ao + aujQ 



bjQO^ + (Ao - po)a - So = 



(35) 



(36) 



with solutions given by a = '^"'^^'^ ■ The exact solutions (1551) . then, follows immediately from (|21l) 

□ 
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Example 2: Consider the following elementary system: 

(j)'^ =301 + 202 (37) 
we have, using ([34| . that a — 1, — |- Thus, for a — 1, the exact solution of the given system is given by 

5 

02 =|Cie4- + C2e-^ (38) 
5 

while for a = — |, we obtain the same solution up to a constant. 
Example 3: Consider the following system 

01 =601 - 02 

0^:^501 +402. (39) 
we have, using (p4)) that a — ^^4^ where i = \/^. The general solution, using ([33]) is given by 



5 

01 ^(i_^)Cie(5+2.). ^ ^^(i_^)^(5-2.). 

02=-^Cie(5+2^)" + C2e(5-2»).. (40) 
which is again agree (up to a constant) with the exact solution obtain by the standard method. 



IV. A CRITERION FOR POLYNOMIAL SOLUTIONS 



In Ref. Saad et al give a sufficient and necessary condition for a second-order linear homogeneous differential 
equation to have a polynomial solution. Although a similar criterion for the existence of a polynomial solution of a 
linear system such as ([3]) is not possible in general, the following theorem gives a necessary condition. 

Theorem 3 If (f)i{x) and 02(2^) are n— degree polynomial solutions of a linear system with variable coefficients 
then 

r]n{x) = Sn{x)uJn{x) - Xn{x)pn{x) ^ 0, n=l,2,... (41) 

where A„,s„,ijj„ and p„ are given recursively by \25]) . 



Proof: The existence of nontrivial solutions 0i and 02 of Q required the vanishing of the determinant 



(42) 



□ 

Note that the converse of this theorem is not true in general, for example, consider the system 0'i = 5a;(0i + 02), 02 = 
3a;(0i +02), clearly rii[x) = 0, however the system has a non-polynomial general solution given by 0i = |Cie'*^ — C2, 
and 02 = fCie^-' + C2. 



V. ELEMENTARY SYSTEMS OF DIFFERENTIAL EQUATIONS 



By means of the iteration sequence (|25|) . it is interesting to note that the termination condition (|22|) can be written, 
equivalently, as Riccati equation, 

{-) ' + (Ao - po) {-)-s, ' = -w, (43) 
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where the solution of this equation yields the exact analytic expression of the ratio a = -^^ by which the exact solution 
of the corresponding system follows immediately using ((2T|) . Thus, we have a useful connection between the exact 
solution of a linear-system of differential equations and the exact solution of a Riccati equation; namely, 



Theorem 4 The general solution of the linear system: 

(f>[ =Ao(a;)(/)i + so(a;)(/)2 

(^2 =UJo{x)(f)i + po (2^)02 

is given by: 



(44) 



4>2ix) = cxp / {po{t) - aujoit)) dt 



=Ci exp 



C2 + C1 / ojoit) exp / (Ao(t) - po{t) + 2ujo{T)a) dr dt 



[aujoit) + Ao(i)) dt - a4>2{x) 



where a = a{x) is the solution of the Riccati equation 



da 2 /X \ 

— ^woa + [Xq - po)a - So. 

ax 



(45) 



(46) 



Clearly, for linear constant-coefficient systems, (1^5)1 is equivalent to (|34p because da/dx = 0. Recently, Saad et al 
studied the exact solutions of certain classes of such Riccati equation which allow us to obtain exact solutions to 
many systems of linear differential equations. We, first, consider the simpler case, namely, the first iteration using 
Theorem 1, we obtain: 



Theorem 5 For the following system of differential equations: 

(\)\ =Ao(x)(/)i +so(a;)(/>2, 
(/)2 =wo(x)(/)i + /5o(a;)(/)2, 

if the functions Aq, Sq; '^O; 01,'f^d po satisfy: 



Ao\' 

so. 



Po 



-LOq 



(47) 



(48) 



then the general solution to the system is given by: 



/,2 

sa{x) 
61 =Ci exp 



wa{t)so{t) 



Xo{t) 



■ exp 



Ao(t) + 



Ao(r) 



dT\dt + C2 



Xoit) 



so{t)wo{t) 



dt 



sojx) ^ 
\q(x) 



Proof: Equation (|^ can be written as 



d 
dx 



Wq 



(Ao - Po) So 



For n — 0, a= j^, and ([50)1 reduces to 
As direct examples of Theorem 5, we have the following two corollaries. 
Corollary 1 The following system of differential equations: 

0'l = f{^) i4>l + <t>2) 
<i'2 = 9{^) (01 + 02) 

where f(x) and g{x) are arbitrary functions, has the general solution given by: 



(49) 

(50) 
□ 



(51) 



62 = Ci 



exp 



{fit)+g{t))dt - 



git) exp { / ifiT)+g{T))dT dt 



-C2, 



5(<)cxp / (/(r)+g(T))dT dt + Q 



(52) 
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TABLE III: Exact solutions of the system (j)'-^ = \o[x)cj>i + so(a;)<^2, ~ uio{x)4>i + po{x)(j}2 for certain functions 

Ao(x), so(2;), uio(a;), and pq{x). 



Xq{x) 


so{x) 


ojq{x) 


po{x) 


Solution 


f(x) 






f(x) 


X _ (7 g/"(9(t)«(t) + /(t))dt _ ^t^\M 

where a(x) = tanh(— g(t)dt). 


fir) 


ag{x) 




f(r\ 
J y-'^i 


Jj^ _ (J S''(.bg(t)c(t) + f(t))dt _ 

^2 = e/"(*)-'«(')"('»''*[C2 + Ci r &5(t)e^'(''"'("'"("»''"dt], 
where a{x) = — ^/|'tanh(\/afe g{t)dt). 


f(x] 






f(x) 


where a{x) = tan(J^ g{t)dt). 


fix) 


-ag{x) 


bg{x) 


fix) 


4>i = Cie-^'('"'(*'"(*'+^('»* - a(:c)<^2, 
where 0(2;) = ^/|'tan(\/a6 g{t)dt). 


Xo{x) 


so{x) 


UJo{x) 


d 1 po-Aq \ 
dx y so J 


^,-C,ei"''^o+^o)dt (Ao-Po)^^, 

[(po- ""'P )dt.„ , „ f'(Ao-po+ 


\o{x) 


d ( Aq-po \ 
dx V WQ ' 


Loo{x) 


poix) 


<^2 = eS ^°*[C2 + Ci /" woeS*^'"'-^°^''^dt]. 



Corollary 2 T/ie following system of differential equations: 

01 = /(X) (01 - 02) 

02 = 5(2:) (01 - 02) (53) 

where f{x) and g{x) are arbitrary functions, has the general solution: 



61 = c 



exp / ifit) - git)) dt]+ gix) exp / (/(r) - ^(t)) dr ] ] dt 



-Co 



h^C^J l^git) exp (fir) - gir)) dr j j dt - C2 (54) 

Using theorem (4), we can now prove the following. 
Theorem 6 The system of differential equations: 

01 =Ao(x)0i + so(a:)02, 

02 =so(a;)0i + Po(a;)02, (55) 

where the functions Ao, sq, and po are assumed to be continuous and differentiable functions, is analytically solvable 
if the quantity (p(x) — Ao(a;))/so(a^) is independent of x. 

Proof: If the quantity ip{x) — Xq{x)) / sq{x) is independent of x, then the Riccati equation becomes a separable 
equation. □ 

Theorem 7 The system of differential equations: 

01 =Ao(x)0i + so(a;)02, 

02 =wo(a;)0i + Po(a;)02, (56) 

where the functions \q, sq, wq and po o,^e assumed to be continuous and differentiable functions, is analytically solvable 
if the quantities WQ{x)/ip{x) — Ao(a;)) and soix) / (pix) — Ao(x)) are independent of x. 
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TABLE IV: Exact expressions of the ratio a = , n — 1,2, . . . for the system 4>'i = Xo{x)4>i + so{x)4>2, 02 = 'ujo{x)4>i + poix)4>2, 
for different Xo{x),so{x),u!o{x) and po(x). The exact solutions of the system then follows by direct substitution of a in (|45ll . 
Here R{x) is an arbitrary differentiable function. 



Xo(x) 



so(x) 



U)o(x) 



4n 



/ R'(x) 
[WW) 



2x 



-i?(5 



2x 



]_ 



R(x) 



R'{: 



R(3 



2(2n + l) 
R{x) 



( R'(^) 
[r^ 



2x 
R{x) 



'R{x 



2x 



]_ 

R{x 



yx a -i\^ — n;^;x" ) -il^j.^ y 

—2n(ax + b] ^- — h n — ax — o 

V ^ ^ iFi(-n;l;lii::#ir /_ 



2na 



/ R'(x- 



I _ (ax + b) 



-i?(a; 



i_ 

R(x 



R.'(x 



R(x 



(2n + l)a 



R{x) 



R(x) 



-R{x 



ax + b 



]_ 

R{x 



ax-\-b 



2n(ax + b) lFl(-n + l;f ; '°"^+''' ) 
3 ,p,^.„.3.,l^^- 



+ 



R(x) 



R(x 



nb 



xR(x) 



R'(x) _ (fc-f ) 
R^x) R(x) 



-i?(5 



b- 



1_ 

R(x 



nb iFi(-n + l;c+l;bx) . R' (x) . 
c iFi(-n\c;bx) ' R(x) ' 
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R'{x 



x(l~x)R(x) 



yww) 



(-2n + l)x- 
x{\-x)R(x) 



-R{x 



( — 2n + l)x — c 
x{l — x) 



1_ 

R(x 



2Fi{-n + l,-n + l;c+l:x) 
2 Fi {—n, — n;c;x) 



+ 



R'(x) _ (-2n+l) 
R(x) 



x — c \ 



R'{x 



R(x 



nb 



x(l-x)R(x) 



R'(x) 
WIT) 



(-n + b+l)x- 



x(l-x)R(x) 



( — n + 6 + 1 ) x — c 
x{l — x) 



]_ 

R(x 



_nb 2-Fi(-n + l,b+l;c+l;3:) 
c 2 -fl i — n,b-c;x) 



+ 



R'(x) 
R(x) 



( — n + b+l)x — c 
x(l — x) 



{ax + b) \' ' 
x'-'Rix) I 



R'{: 



R(3 



n(n + a — 1) 
x'^R{x) 



R'jx) 
WIT) 



-R{i 



ax-\-b 



1 

'Wx) 



n{n + a-l) 2 Fp ( -n+l,n + a; - 

b 2-F'o(-"."+a-l;- 



ll I -R'(3;) (gjj+b) 



R'[x 



R(3 



n(n + 2k + l) 
(l-x'^)R(x) 



R'jx) 
WIT) 



2(k + l)x 
(l-x'-')R{x) 



-R{x 



2(k + l)x 



1_ 

R{x 



n(n+2fc + l) 2Fi(-n+l,n+2fc + 2;fc + 2;ii-^) _^ ^'(j:) 



2{fc + l) 



2Fi{-n,n + 2fc+l;fc + l;-^ 



R(x) 



2(k + l)x 



n(u+2fc) 



Rf(x)_ 
WIT) 



(2k + l)x 
l-x''')R{x) 



-i?(5 



(2fe+l)£ 



]_ 



n(n+2fc) 2Fi(-n + l,n + 2fc + l;fc+f 



(2fc + l) 



2Fi{-n,n+2k;k+} 



(2fe + l)£ 



Proof: If the quantities wo{x) / {p{x) ~ Xo{x)) and sq{x) / {p{x) — Xq{x)) are independent of then the Riccati equation 
becomes a separable equation. □ 

VI. SOLVABLE LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS 

For higher iteration levels, n — 1,2, ... , the results of the earlier work of Saad et all @ on Riccati equation [T^ - 
[3] can be used to obtain exact analytical solutions to different classes of the linear system (|3]). Under a certain 
conditions on the functions Aq, sq, wq and po, we may express [l5j | some of these exact solutions in terms of generalized 
hypergeometric functions qFp{ai, . . . , ap; Pi, . . . , I3q;x) 



pFq{ai,a2, . . ., Up] Pi, (32, . . . , a;) = ^ 



(ai)fc(a2)fc ■ ■ ■ (ap)fc x^ 

{Pl)k{f}2)k...{Pq)k k\ 



(57) 



fc=0 



where p and q are nonnegative integers and no A; = 1, 2, . . . , g is zero or a negative integer. Clearly, (I57|) includes 
the special cases of the confluent hypergeometric function iFi and the classical 'Gaussian' hypergeometric function 
2F1. The Pochhammcr symbol (a)fc is defined in terms of Gamma function as 



, , Via + n) 

r(a) 



(58) 



If a is a negative integer —n, we have 



r (-1)""! < A- < 77 
(-n)fc = ("-'=)' ^ ^ ^ « 
k> n 



(59) 



in which case, the generalized hypergeometric series reduces to a polynomial of degree tt, in its variable x. In Table 
IV, we report the exact expression of a for a given system, by which we can easily obtain the exact solutions through 
a direct substitution of a in (1211). 



Theorem 8 If functions Fq = Fq{x) and Gq = Gq{x) satisfy the recursive relation: 



FnGn-i - Fn-iGn = 0, for somc 



77= 1,2,. 



(60) 
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where 

\Gn = G[^_i + Goi^n-1, 

then the linear system of differential equations: 

=Ao(/)i + Go exp ^- J {Fo + po - Ao)dt^ 02 
0'2 = exp ( / (Fo + po - Ao)dt ) -/-i + po'/'2 (62) 



has a general solution given by ^5\ ) with Sq — Go exp (— {Fo + Po ^ ^a)dt) , = exp (^J^ {Fq + Po — ^o)dt^ , o,rid 

Gn-l 



^ ^ exp y {Fo + po- Xo)dxj , (63) 
where Xq = Ao(x) and sq = sq{x) are arbitrary functions. 

Proof: The substitution a = — transfer the Riccati equation ((46|) into the second-order differential equation 

u' + sqWqu (64) 
Using the standard theorem of the asymptotic iteration method 3], the solution of satisfies @ 



— + Ao - Po 
Wo 



G„_i 

— = -T^ (65) 

if, for some n > 0, 

5n = FnGn-1 — Fn-lGn = 0, 

where Fn and G„ are computed using the recursive relations 

{Pn = Fn-1 + Gn-1 + FoFn-1, 
Gn = GJj_]^ + GoFn-l- 

Thus, the general solution of (|68p . using Theorem 4, is given by (I45p with 

Wotn-l 

which complete the prove of the theorem. CH 

Example 4: Let Fo — 2x and Go = —2n, n = 1,2 ... . Using (pTjl . see also Q, we know that (5„ = FnGn-i — 
Fn-iGn = for n = 1, 2, . . . with 

G„_i _ 1 dHn{x) 



Fn-l Hn{x) dx 



(67) 



where Hn{x) is the well-known Hermite polynomials. Thus, for po = po{x) and Ao = Ao(a;) arbitrary functions, the 
system 



01 =Ao0i - 2n exp y (2a:: + po - Ao)dtJ 02 

02 = + Po - Ao)dt^ 01 -I- po02 (68) 
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TABLE V; Exact expressions of the ratio a 



^o{^)F„-t , n- - 1, 2, . 



for the system = \o{x)(f>-i + so{x)(f>2, 02 = 



wo{x)(t>i + po{x)<f>2, for arbitrary functions \o{x) and pq{x). The exact solutions of the system then follows by direct substitution 
of a in (l45ll. 



so{x) 



LOo{x) 



1 G,^^l 



-2ne' 



-/="{2a: + po-Ao)da: 



,/^{2a: + po-Ao)da; 



J"^ (ax + b+pQ-Xa)dx 



^J"^ (ax + b+pa-\Q)dx 



lFi(- 



n 1 (aic + by_ 
Jl±22 



- {ax + h+p^i — XQ)dx 



f"^ (ax + b+pQ-Xa)dx 



n = 1, 3, 5, 



,f'^ (ax + b+pQ-\Q)dx 



„-1.3, (ax + 6)^ 
2 ' 2 ' 2a 



{ax + b)iFi{- 



{ax + b+ PQ — \Q)dx 



_bri^- {b- %+Po->'o)<lx 



J''{b-i+PO->'o)<lx 



1-^1 ( — 'n';c;63;) 



j^bFl(-n,-n;c;x)] ^-J=^(ill^ 7+_l),^-'- +P0-A;r^ 



— x) 



j^[2Fl(-n,n + l;l;ii^)] - ( +po - Ap ) 



n{n + l) ■ 



-+P0-A0 da; 



r+PO--'>0 da; 



lFi(-n,Ti + l;l; ''-^°'' ) 



n(n + a + t.+ l) „ 5^ ( 



{l-a;2) 



(a + fc + 2)x + a-b ^ 

TT^;:^^ +P0-A0 , 



(a + 6 + 2)x + a-6 



j^[2Fl(-"," + a + ft+l;a + l;-^^)J -/ 
2Fi(-n,Ti + a + 6+l;a + l; ) 



(a + b+2)x + a-b 



+ P0-^0 da; 



n(n + 2) „- J^(tZ^+PO~^o)° 



J"(T^+PO-Ao)da; 



2Fi( — n,n+2; 



2J 2_ 



n(n + 2fc) ^-J"(^^;±^+PO-Ao)da; 



-f!- 2Fi(-n,n+2fc;fc+J 



2Fi(-n,n + 2fc;fc+^;ii^50 



jL^2Fi(-n,n + 2fc + l;fc + l;ii-^)J ^_ ( +p„ _ Aq) „ 

2Fi(-ra,n+2fc + l;fc + l; ^^1,"^ ) 



(n+2fc + l) - J°'( '''/l''J^2°' +PO-^o)° 



ji^+Po-Ao)d. 



J^-^+Po-^^ 



n(n + l) ^- J" (^Hiij^+po-Aoy: 



s[2f'o(-"."+i;--f)j -r(-^^+PO-Ao). 



n(n + a-l) ^- J" ( '7/"' +PO-Ao)c 



J"(-i^+Po-Ao)da; 



2Fo(-n,n+l;-;-§) 



E[2fo(-"."+"-i;-;-f)j ^-r(-^^+Po-^ 



2 Fq { — n,n-\-a — \; 



has a general solution 



'^i(^) ^irfr^^ "°^'^'''-a02(x). (69) 

G 

In Table V, we gave the exact expressions of a = ^^^J^"'^ ^ , n = 1, 2, . . . for different systems. The general solutions 



can be found by direct substitution of a in (j45|. 

Theorem 9 If functions Fq = Fq{x) and Gq = Gq{x) satisfy the recursive relation: 

Sn = FnGn-i - F^-iGn = 0, for souic n = 1, 2, . . . (70) 

where 

{Fn = F'^_i + Gn-l + -Fo-F'n-l, ^rj^-^ 
Gn — + 

then the linear system of differential equations: 

4>'i =Ao0i + exp (Fq - po + Ao)da;^ 02 

02=G'oexp^y" -(-Fo - Po + Ao)rfa;^ 01 + Po02 (72) 
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TABLE VI: Exact expressions of the ratio a ; 



soF„-l 
G„-i 



, n = 1,2, . . . for the system = \o{x)cj>\ + ao{x)cj>2, 02 = wo(a;)</)i + 



po{x)(t>2, for arbitrary functions Ao(a::) and po{x). The exact solutions of the system then follows by direct substitution of a in 
(|45]l. 



So(a::) 



1^0 (a;) 



-2ne~ ■'^^ (23;-po+Ao)da: 



(aa; + 6-po + Ao)tia; 



n = 2,4,6, 



lJ"l(-2-'2- 2a ) „f"(a3: + 6-po + Ao)da: 



/^(ai + 6-po + -'io)da; ^ _ g g 



(a3; + i))lFl(-:g^;f ; <°°^+^''' ) (ax + b- pa + \a)dx 



i (a. + 6)iFi(-^;|;<2^) 



lFl(-n;c;63:) f ^ ( t- ^ -pp +An ) da: 
_£[l£l(-n^5Ml 



j|[2Fl(-n,-n;c;3;)] 



x{l — x) 



J'' ( (i^^2 -~PO+>-oj~ 



lFi(-n,n+l;l;i 



2Fi(-n,n+l;l; 



t1 , 1 



rx I (a + b + 2)x + a-b , \ , 

(^^^ (i_!,2) PO + ^0 J da: 



rx / (a + 6 + 2)x + a-b „ , , \ 



2Fi(-n,n + a + b+l;a + \;^-^-^) J'' 



(a + 6 + 2)x + a-6 



-po + Ao da: 



(YTfT-PO+Ao)da: 



n(n + 2) „- /"(tTTT-PO + ^o)''^ 



2Fi(-Ti,n+a + h+l;a + l; ) 



.3. (1-x) 



2Fi( — n,n+2; 



,r(^-po+Ao)d.- 



(n+2fc) -J-(i^^-po + Ao)d.- 



^[2Fi(-n,n+2;i;-^i^ 



-PO + Ao Ida: 



2Fi(-»,n+2fc;fc+^;il^) J" ( ^^^^^ -pp +Ao ) da: 
2Fi(-n,n + 2fc;fc+i;iij^) 



.r(^¥^-po+A.;)^ 



n(„ + 2fc + l) - J"(^il-^-po + Ao)da:~ 



2Fi(-n,n + 2fc + l;fc + l;i^-^) J ( ^ '/l^j'j -po + Ao)dj: 



(l-a:2) 



,r(^^-po+A";;)d. 



^ 2f'i(-n,"+2fc+l;fc+l;^^-2^) 



2Fo(-n,n + l;-;-f ) ^/^ ( - ^^1 + ") -pp + Ap ) da: 
£[2Fo(-n,u+l;-;-f)] 



.J"( ''r2"^'" ~PO + Ao)da: 



n(n + a-l) ^- J°' ( - i^l+ii -po + Ap )~ 



2Fp(-n,n + a-l;-;-f) ( - i£i£j±i -pp +Ap ) da: 
j^[2Fo(-n,n + a-l;-;-f)] 



/las a general solution given by wii/i so = exp [J^{Fq — po + Xo)dx) , wq — exp (— J^{Fq — po + Xo)dx) , and 



Fn- 



- exp ( {Fo + Pa ~ Xo)dx] , (73) 



Gn-l 

where Xq = Xq(x) and sq = so{x) are arbitrary functions. 

Proof: The proof of this theorem follows similarly to the prove of Theorem 6 using the substitution a = — ^ . D 

In Table VI, we gave the exact expressions of a = ^jr-^^, n = 1, 2, . . . for different systems. The general solutions 
can be found by direct substitution of a in ()45|) . 

VII. CONCLUSION 

Guided by the recent applications of the asymptotic iteration method on Dirac equation and Riccati equation [3]-[6[, 
we extend AIM to solve analytically different systems of linear differential equations with variable coefficients. This 
allow us to provide a simple recipe for solving any linear constant-coefficient systems in few steps. Further, we have 
extended the standard list of exactly solvable system known in the literature (see [ij, Chapter T6, Section T6.1.1, Page 
1229). Furthermore, the work presented here established a connection between the solution of arbitrary linear system 
of differential equation and the solution of Riccati equation. Indeed, one of the main results in the present work, is 
to have this connection written explicitly in terms of a theorem; namely Theorem 4. Another interesting conclusion 
of the present work is that it provides a set of conditions on Xq, so,wo, po which govern the complete solvability of 
the system ([3]). It should be clear the results of our investigation are not limited to the exact solutions reported in 
tables I- VI but it can easily extended to many other solvable systems. General speaking, the results obtained here 
are useful in the sense that they can set up a more rational base to approach to solve coupled differential equations 
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[l6j-|17l|. provide a means to examine this problem from a qualitative point of view, alleviating then the burden of 
the computational work. 
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